home *** CD-ROM | disk | FTP | other *** search
/ Brotikasten / BROTCD01.iso / lern / lern338.d64 / ma010 bairs eq (.txt) < prev    next >
Commodore BASIC  |  1995-08-09  |  1KB  |  67 lines
  1. 1 REM *** BAIRSTOW
  2. 2 REM     W HAWES
  3. 3 REM
  4. 4 REM *** PET USER GROUP
  5. 5 REM     SOFTWARE EXCHANGE
  6. 6 REM     PO BOX 371
  7. 7 REM     MOTGOMERYVILLE,PA 1936
  8. 8 REM
  9. 10 POKE59468,12
  10. 19000 PRINT"[147]"
  11. 19010 REM PROGRAM SUBMITTED BY W. HAWES
  12. 19020 REM HALIFAX, N.S. CANADA
  13. 19030 PRINT"PROGRAM BAIRSTOW"
  14. 19032 PRINT
  15. 20000 PRINT"PROGRAM TO SOLVE AN N'TH ORDER"
  16. 20002 PRINT"POLYNOMIAL EQUATION OF THE FORM"
  17. 20004 PRINT"Y^N + A1*Y^(N-1) +....+AN = 0
  18. 20006 [153]"WHERE A1,A2,....,AN ARE REAL.
  19. 20008 PRINT"BAIRSTOW'S ITERATIVE METHOD IS
  20. 20010 [153]"USED TO YIELD SUCCESSIVE
  21. 20012 PRINT"QUADRATIC FACTORS AND THE
  22. 20014 [153]"RESULTANT REDUCED ORDER
  23. 20016 PRINT"POLYNOMIAL
  24. 20020 [153]
  25. 20022 [153]"PRESS ANY KEY WHEN READY"
  26. 20024 [161]A$:[139]A$[178]""[167]20024
  27. 20100 [153]"LOAD"
  28. 20102 [133]"ENTER ORDER N OF POLYNOMIAL";N
  29. 20110 P[178]0:Q[178]0:Z[178]0
  30. 20120 [134]A(N),B(N),C(N):C(0)[178]1
  31. 20130 [129]L[178]1[164]N
  32. 20140 [153]"ENTER A";L;:[133]A(L)
  33. 20150 [130]L
  34. 20155 [153]:[153]
  35. 20160 B(1)[178]A(1)[171]P:B(2)[178]A(2)[171]P[172]B(1)[171]Q
  36. 20170 [129]K[178]3[164]N
  37. 20180 B(K)[178]A(K)[171]P[172]B(K[171]1)[171]Q[172]B(K[171]2)
  38. 20190 [130]K
  39. 20200 C(1)[178]B(1)[171]P:C(2)[178]B(2)[171]P[172]C(1)[171]Q
  40. 20210 [129]K[178]3[164]N
  41. 20215 C(K)[178]B(K)[171]P[172]C(K[171]1)[171]Q[172]C(K[171]2)
  42. 20220 [130]K
  43. 20230 CB[178]C(N[171]1)[171]B(N[171]1)
  44. 20240 DP[178]B(N[171]1)[172]C(N[171]2)[171]B(N)[172]C(N[171]3)
  45. 20245 DP[178]DP[173](C(N[171]2)[172]C(N[171]2)[171]CB[172]C(N[171]3))
  46. 20250 DQ[178]B(N)[172]C(N[171]2)[171]B(N[171]1)[172]CB
  47. 20255 DQ[178]DQ[173](C(N[171]2)[172]C(N[171]2)[171]CB[172]C(N[171]3))
  48. 20260 P[178]P[170]DP:Q[178]Q[170]DQ
  49. 20270 M[178][182](DP)[170][182](DQ):Z[178]Z[170]1
  50. 20280 [153]"ON";"ITERATION #";Z;" ERROR =";M
  51. 20290 [139]M[179].0000001[167]20320
  52. 20300 [137]20160
  53. 20320 [153]
  54. 20340 [153]"P = ";P
  55. 20350 [153]"Q = ";Q
  56. 20360 [153]
  57. 20370 [153]"P & Q ARE A QUADRATIC FACTOR OF"
  58. 20380 [153]"THE FORM Y^2 +P*Y + Q = 0"
  59. 20390 [153]
  60. 20400 [153]"THE RESULTANT FACTOR IS OF THE"
  61. 20410 [153]"FORM: Y^(N-2) + B1Y(N-3)+...+B(N-2)"
  62. 20420 [153]"WHERE:"
  63. 20430 [129]K[178]1[164]N[171]2
  64. 20440 [153]"B";K;" = ";B(K)
  65. 20450 [130]K
  66. 20460 [128]
  67.